Effects of the field modulation on the Hofstadter's spectrum 
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We study the effect of spatially modulated magnetic fields on the energy spectrum of a two- 
dimensional (2D) Bloch electron. Taking into account four kinds of modulated fields and using 
the method of direct diagonalization of the Hamiltonian matrix, we calculate energy spectra with 
varying system parameters (i.e., the kind of the modulation, the relative strength of the modulated 
field to the uniform background field, and the period of the modulation) to elucidate that the energy 
band structure sensitively depends on such parameters: Inclusion of spatially modulated fields into 
a uniform field leads occurrence of gap opening, gap closing, band crossing, and band broadening, 
resulting distinctive energy band structure from the Hofstadter's spectrum. We also discuss the 
effect of the field modulation on the symmetries appeared in the Hofstadter's spectrum in detail. 

PACS numbers: 71.28.+d, 71.25.-s, 73.20.Dx, 71.45. Gm 



I. INTRODUCTION 

The problem of a 2D Bloch electron under a uniform 
magnetic field has been intensively studied for several 
decades, 1 and it is well known that the energy spectrum 
is characterized by the Hofstadter's butterfly showing a 
fractal nature. 2 Recently the problem has attracted theo- 
retically renewed interest in connection with various phe- 
nomena such as the quantum Hall effect, 3 the flux state 
model for high-T c superconductivity, 4 and the mean-held 
transition temperature of superconducting networks or 
Josephson junction arrays. 5 Besides, recent advance in 
submicron technology that makes it possible to fabricate 
any desired microstructures has stimulated experimen- 
tal study to find indication of the Hofstadter's spectrum 
and its effect on the transport and optical properties. 6-9 
In parallel with this problem, the study on 2D elec- 
tron systems under nonuniform (either disordered 10 or 
periodic 11-26 ) magnetic fields has been extensively per- 
formed, and lots of interesting characteristics in the en- 
ergy spectral and transport properties have been eluci- 
dated. 

Though the problem of a 2D Bloch electron under spa- 
tially modulated magnetic fields has attracted less atten- 
tion compared with the problem of 2D electron gas un- 
der spatially modulated magnetic fields, 11-22 it is still 
an important problem not only in the view of theoreti- 
cal interest but also in the view of experimental interest, 
and there have been attempts to solve this problem. 23-26 
However, unfortunately, some of the relevant works con- 
tain inconsistent results on the energy spectral proper- 
ties: In Ref. 24, Gumbs et al. studied the effect of a 
one-dimensional sine-modulated (1DSM) field on a 2D 
Bloch electron to argue that the symmetries appeared in 
the Hofstadter's spectrum break down by the field mod- 
ulation. Most surprisingly, they also argued that the 
field modulation leads an additional crisscross pattern 
like a spiderweb structure onto the Hofstadter's spec- 
trum. However, Oh et al. 25 studied the effect of a ID 



cosine- modulated (1DCM) field on a 2D Bloch electron 
to elucidate that occurrence of gap closing and gap open- 
ing leads different energy spectrum from the Hofstadter's 
spectrum and that the symmetries appeared in the Hof- 
stadter's spectrum except the dual property still remain 
despite the field modulation. In the meanwhile, Shi and 
Szeto 26 studied the energy spectrum of a 2D Bloch elec- 
tron under a kind of 2D field modulation to argue that 
there is no symmetry breaking in the energy spectrum 
and that the fractal structure remains irrespective of the 
field modulation. 

In this paper we reexamine the problem of a 2D Bloch 
electron under spatially modulated magnetic fields to set- 
tle the inconsistency discussed above. In doing this, we 
take into account four kinds of modulated fields (i.e., 
1DSM, 1DCM, 2DSM, and 2DCM fields) in order to ob- 
tain rather generic effects of the field modulation on the 
Hofstadter's spectrum. By means of direct diagonaliza- 
tion of the Hamiltonian matrix, we calculate the energy 
eigenvalues and examine how the system parameters such 
as the kind of the modulation, the relative strength of the 
modulated field to the uniform field, and the period of 
the modulation influence on the energy band structure 
and the symmetry of the Hofstadter's spectrum. Intro- 
duction of the field modulation is shown to change the 
k dependence of the energy spectrum drastically, leading 
occurrence of gap opening, gap closing, band crossing, 
and band broadening, which are the origin of distinctive 
energy band structure from the Hofstadter's spectrum. 
Our results indicate that there is no additional spider- 
web structure in the energy spectrum contrary to the 
result of Ref. 24 and that the field modulation generi- 
cally breaks the symmetries and the fractal property of 
the Hofstadter's spectrum. 

This paper is organized as follows: In Sec. II we intro- 
duce four kinds of magnetic fields and the tight-binding 
model. And, in Sec. Ill we present numerical results on 
the energy band structure of a Bloch electron and the 
symmetries of the energy spectrum under these fields. 
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Section IV is devoted to a brief summary. 



II. MODULATED MAGNETIC FIELDS AND THE 
TIGHT-BINDING MODEL 
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We consider an electron in a 2D square lattice under a 
spatially modulated magnetic field 



B = [B + B 1 (x,y)]z, 



(1) 



where Bg (Bi) denotes the uniform (modulated) part of 
an applied magnetic field. Among possible kinds of mod- 
ulated fields, we pay attention to two kinds of modulated 
fields. The one is the SM field 



B^x.y) = B x sm 



/2nx\ 
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and the other is the CM field 

' 2nx 



B^x.y) = B x cos 
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Here B x / y \ is the strength of the modulated field and 
T x fy\ is the period of the modulation along the x (y) di- 
rection. Under the Landau gauge, the vector potential 
becomes 



ByTy 
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for the SM field and 



\T V J' 



Ay = B x 
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for the CM field, respectively. 

The tight-binding Hamiltonian describing the motion 
of an electron in a magnetic field is given by 



(6) 



where tij is the hopping integral between the nearest- 
neighboring sites i and j, and dij = (2ne/hc) J. A ■ dl 
is the magnetic phase factor. Under the vector poten- 
tials given by Eqs. (^) and (||), the magnetic phase factor 
becomes 

±9 m , (m , n ) = (m, nil) 
±6 n , (m ,n ) = (m± l,n) (7) 
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for the SM field and 
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for the CM field, respectively. Here fi x ( y ) = B x ^/Bo, 
lx(y) = T x ( y )/a, and <po = Bqo 2 , a being the lattice con- 
stant. The magnetic flux per unit cell, in units of the flux 
quantum hc/e, is given by = (l/27r) ^2 &ij ~ § A ■ dl — 
jB-dS. 

By means of Eqs. (||) and (]?]), the tight-binding equa- 
tion can be written as 



-i6„ 



4>r, 



(10) 



where A (= t y /t x ) is the ratio of hopping integrals be- 
tween the x and y directions, and E is the energy in units 
of t x . Here, the wave function is given by \ip) = • ipj\j)- 

Denoting R x ( y ) as the periodicity of 9 m (n), the Bloch 
theorem can be written as 



Ipmn, 1pm,n+Ry = < 5 ' ' V 1pmn, (H) 



and the first magnetic Brillouin zone (FMBZ) is given by 
%/R x / y y We calculate the energy eigenvalues 



\k. 



< 



for all the values of k in the FMBZ by directly diago- 
nalizing the Hamiltonian matrix obtained from Eqs. (|10| ) 
and (O). 



III. NUMERICAL RESULTS AND DISCUSSION 



In what follows, we assume the modulated field has a 
square-symmetry (i.e., (3 X — (3 y — (3 and j x = j y = 7) 
and consider only the case of the isotropic hopping in- 
tegral (i.e., A = 1) for the sake of simplicity. And we 
pay attention to the energy dispersions for q = 2, 3 (with 
setting p = 1) and 7 = 2, 3, 4 because energy dispersions 
for other values of (q, 7) can be obtained in a similar way. 
Here p and q denote the numbers (prime each other) given 

by 0o = p/q- 
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A. Energy spectrum in a uniform magnetic field 

In a uniform field ({3 = 0), 9 n becomes zero. Thus, by 
means of the translational invariance along the y direc- 
tion, Eq. (|l(i| ) reduces to the Harper equation 



ipm+i + ipm-i + 2cos(6» m + k y )ip m = Eip T: 



(12) 



and there are q eigenvalues for a given value of k y . When 
q is odd, the full energy spectrum consists of q subbands. 
In the meanwhile, when q is even, the full energy spec- 
trum consists of (q — 1) subbands because two central 
subbands touch at zero energy. Here the full energy spec- 
trum means a set of energy eigenvalues obtained by tak- 
ing into account all the values of k in the FMBZ. Energy 
dispersions for q — 2 and 3 are given as follows: 

(a) q = 2: We have R x = 2 and the energy dispersion 
is given by 



E(k) = ±2Wcos 2 k x + cos 2 k y 



(13) 



There are two subbands and they touch at zero energy as 
shown in Fig. 1(a), resulting a single subband structure 
with \E\ < 2V2. 

(b) q — 3: We have R x = 3 and the energy dispersion 
is given by the equation 



E 3 - 6E-2(cos3k x 



cos3k y ) = 0. 



(14) 



There are three subbands as shown in Fig. 1(b). But, 
there is no touching point in the dispersion and the 
full energy spectrum exhibits a three subband struc- 
ture with -(1 + VS) < E < -2, \E\ < (V3 - 1), and 

2 < E < (1 + \/3)- 

It is well known that the Hofstadter's spectrum has 
the following symmetries; (i) the dual property between 
the k x and k y directions (let us denote it as Sd), (ii) 
the symmetry between —E and E (Se), (hi) the sym- 
metry between — k x and k x (Sx), and (iv) the symmetry 
between — k y and k y (Sy)- It should be noted that Se 
holds only when all the values of k in the FMBZ are taken 
into account; for a given value of fc, one can easily check 
that Se does not hold for odd q. 



B. Energy spectrum in ID modulated magnetic 
fields 

When the held modulation is along the x direction, n 
is still zero and the tight-binding equation is formally the 
same as Eq. (12). However, [3 and 7 are introduced in 
m , and (j) becomes periodic along the x direction with 
the period 7: 



1 + (f3j/2n) [cos (2tt(to- 
cos (27rm/7)] 



i)h) 

for 1DSM 



1 - (/37/27I-) [sin (27r(m + IV7) 

- sin(27rm/7)] for 1DCM. 



(15) 



The lattice is called the stripped flux lattice 23 and the 
energy spectrum can be obtained by diagonalizing the 
(Rx x R x ) Hamiltonian matrix. 



1. 1DSM Field 

When (9,7) = (2,2), we have R x 
dispersion is given by 



2 and the energy 



E(k) = 2 sin (3 sin k y ± 2 y cos 2 k x + cos 2 (3 cos 2 k y . (16) 

The two subbands touch at (k x ,k y ) = (ir/2,ir/2) and 
(ir/2, 3ir/2), and the energy eigenvalues at the touching 
points are given by E — ±2 sin j3. Since the touching 
points exist irrespective of (3, the full energy spectrum ex- 
hibits a single subband structure as in the case of [3 = 0. 

In order to demonstrate the 7 dependence of the en- 
ergy spectrum with q — 2, we calculate energy spectra for 
various values of 7 and plot some of them in Fig. 2. Even 
though there are six and four subbands for (q, 7) = (2, 3) 
and (2, 4), there occurs direct touching between the near- 
est neighboring subbands, which leads a single subband 
structure of the full energy spectrum. In the calculations, 
we hnd that the energy band structure is independent of 
the values of (3 and 7 while the total bandwidth slightly 
changes with varying (3. 

Figure 3 shows the k y dependence of the energy spec- 
trum for (q,"f) = (3,2). In this case, we have R x = 6 
and there are six subbands. For small values of (3, the 
upper (lower) two subbands touch at several points of k y . 
Besides, there exists indirect overlapping (i.e., crossing of 
subbands at different values of k y ) between the central 
subbands. Thus, the full energy spectrum still exhibits a 
three subband structure as can be seen in Fig. 3(a). How- 
ever, as (3 increases, there occurs gap opening between 
the subbands and the full energy spectrum exhibits a six 
subband structure as in Fig. 3(b). Further increase of (3 
makes the second and the fourth gaps closed, leading a 
four subband structure as in Fig. 3(c). Figure 4 shows the 
k y dependence of the energy spectrum for (9,7) = (3,3). 
In this case, we have R x — 3 and there are three sub- 
bands. For small values of /3, the full energy spectrum 
exhibits a three subband structure as in Fig. 4(a). How- 
ever, indirect overlapping between subbands occurs as (3 
increases, and the full energy spectrum exhibits a single 
subband structure as in Fig. 4(b). Occurrence of gap 
closing due to indirect overlapping is also found in the 
case of (q, 7) = (3, 4). 

As for the symmetry of the energy spectrum, Figs. 2 
and 3 show that Se remains under the 1DSM held, which 
is contrary to the arguments of Ref. 24. The reason of 
this discrepancy lies in the range of k taken into account 
in discussing the symmetry of the energy spectrum. In 
Ref. 24 only a particular value of k was taken into ac- 
count, while all the values of k in the FMBZ are taken 
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into account in this paper. Here we would like to stress 
that all the values of k in the FMBZ should be considered 
in order to discuss the symmetry of the energy spectrum. 
This is because Se of the Hofstadter's spectrum holds 
only when the all the values of k in the FMBZ is consid- 
ered. We also find that Sx still remains under the 1DSM 
field while whether Sy remains or not depends crucially 
on q and 7. Note that Sd breaks down by introducing 
the ID field modulation. 



2. 1DCM Field 

In order to test how the energy spectrum is influenced 
on the kind of the field modulation, we also calculate 
the energy spectrum of a 2D Bloch electron under the 
1DCM field. Figure 5 shows the k y dependence of the 
energy spectrum for (q, 7) = (2, 3). In this case, we have 
R x — 6 and there are six subbands. For small values of (3, 
the two central subbands directly touch at two points of 
k y and the upper (lower) two subbands indirectly overlap 
with each other, resulting a single subband structure as in 
Fig. 5(a). However, as /3 increases, indirect overlapping 
between the upper (lower) two subbands disappears and 
there occurs gap opening between them. Thus the full 
energy spectrum exhibits a three subband structure as in 
Fig. 5(b). Further increase of j3 makes the remaining indi- 
rect overlapping disappear and the full energy spectrum 
exhibits a five subband structure as in Fig. 5(c). The (3 
dependence of the energy spectrum for (9,7) = (2,4) is 
quite different from (5,7) = (2,3). In this case, we have 
R x = 4 and there are four subbands. Each subband di- 
rectly touches with its neighboring subbands and the full 
energy spectrum exhibits a single subband structure re- 
gardless of j3; the only effect of f3 is to change the energy 
bandwidth. 

We also perform similar calculations for the energy 
spectrum with q = 3. When (9,7) = (3,3), we have 
R x = 3 and there are three subbands. For small values 
of (3, the three subbands separate each other, exhibiting a 
three subband structure as in the case of j3 = 0. However, 
as (3 increases, each subband overlaps indirectly with its 
neighboring subbands and the gaps between the sub- 
bands become closed, resulting a single subband struc- 
ture. We also find similar phenomena for (9,7) = (3,4). 

In Table 1, we summarize the number of subbands of 
the full energy spectra for the parameters (q, 7,/3) we 
took into account. As for the symmetry of the energy 
spectrum, we find that Se, Sx, Sy remain, while So 
breaks down under the 1DCM field. Note that Sy de- 
pends on the kind of the field modulation; it remains 
(breaks) in the 1DCM (1DSM) field. 

Before ending this subsection, we would like to men- 
tion two comments. First, we assume 7 to be an integer 
for the sake of simplicity even though it can be arbi- 
trary real value. Besides, we consider only the values of 
7 > 3 (2) in the case of the CM (SM) field because <t> be- 



comes nonuniform only under these conditions. However, 
the authors of Ref. 24 violated this condition; they chose 
7 = 1 in their calculations. But, when 7 = 1, </> becomes 
uniform as can be easily checked by Eq. ( |l5| ) and the en- 
ergy spectrum should be identical with the Hofstadter's 
spectrum if all the values of k in the FMBZ are taken into 
account, which implies that the results of the cited pa- 
per, i.e., the appearance of additional crisscross pattern 
like a spiderweb structure and the symmetry breaking 
of the energy spectrum, are erroneous ones coming from 
mistakes in choosing the values of 7 and in choosing the 
values of k. Note that the method of diagonalization of 
the Hamiltonian matrix used in the present paper or in 
the cited paper and the method of the transfer matrix 
used in Ref. 25 are equivalent. Actually, Fig. 5(a) in the 
present paper is exactly the same as Fig. 2(b) in Ref. 25. 
Second, the tight-binding model we are considering is ba- 
sically a one-band model. Thus, we focus our attention 
on energy spectra only for the values of (3 that are not 
large (i.e., (3 < 1) since there might be interband mixing 
between different Landau levels for large values of (3. 



C. Energy spectrum in 2D modulated magnetic 
fields 



When the field modulation is along both lateral di- 
rections, 9 n becomes nonzero and <f> becomes periodic 
in both lateral directions with the period 7. The lat- 
tice is called the checkerboard flux lattice 23 and the 
energy spectrum can be obtained by diagonalizing the 
(R x R y x R x R y ) Hamiltonian matrix. 



1. 2DSM field 

When (9,7) = (2,2), we have (R x , R y ) 
the Hamiltonian matrix becomes 



( a b 

a* c 

b* d 

V c* d* 



(2,2) and 



(17) 



where a = -e zl3 - e -W+ 2k y) , b = e" i/3 + e «(^- 2 ^), 

C = e lf3 + e - l (/ 3 + 2fc x) j an( J d — e -ifi _|_ gi(/3-2fc„) j reS p ec _ 

tively. By diagonalizing Eq.(^) we obtain the energy 
dispersion as 



E(k) = ±2 



cos /3\/ cos 2 k x + cos 2 k y 



isin/3-Wsin k x + sin k y 



(18) 



which is plotted in Fig. 6. Numbering the four subbands 
in order of lowering energy, one can see that the first 
(second) and the third (fourth) subbands directly touch 
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at the center of FMBZ regardless of (3. There also ex- 
ists band crossing between the second and the third sub- 
bands, which is absent in the case of the ID field modu- 
lation. Equation (18) indicates that the full energy spec- 
trum exhibits a single subband structure regardless of (3 
and the total bandwidth changes with varying (3. 

When (q, 7) = (2,3), we have (R x ,R y ) = (6,3) and 
there are eighteen subbands. However, due to occur- 
rence of direct touching and indirect overlapping be- 
tween neighboring subbands, the full energy spectrum ex- 
hibits a single subband structure regardless of (3. When 
(9)7) = (2,4), the situation becomes different. In this 
case, we have (R x ,R y ) — (4,4) and there are sixteen 
subbands. For small values of (3, due to direct touching 
between neighboring subbands, the full energy spectrum 
exhibits a single subband structure. However, gap open- 
ing between subbands occurs with increasing (3, resulting 
a five subband structure of the full energy spectrum. 

When (q, 7) = (3,2), we have (R x ,R y ) = (6,2) and 
there are twelve subbands. In this case, even for small 
values of (3, the full energy spectrum is quite different 
from the case of (3 — 0; the full energy spectrum exhibits 
a eight subband structure as can be seen in Fig. 7(a). 
And, as (3 increases, two more gaps are open to yield a 
ten subband structure as in Fig. 7(b). Further increase 
of (3 makes the full energy spectrum exhibit a twelve sub- 
band structure as in Fig. 7(c). When (q,j) = (3,3), the 
effect of (3 is quite different from the case of (q, 7) = (3, 2). 
In this case, we have (R x ,R y ) = (3, 3) and there are nine 
subbands. For small values of (3, the full energy spectrum 
consists of three subbands as in the case of (3 = 0. How- 
ever, there occurs gap closing with increasing (3 and the 
energy spectrum exhibits a single subband structure. We 
also find similar gap closing behavior for (9,7) = (3,4). 



2. 2D CM Field 

In order to know how the energy spectrum is influ- 
enced by the kind of the 2D field modulation, we also 
calculate the energy spectra of a Bloch electron under 
the 2DCM field with various values of (q,j,(3). And, we 
find that occurrence of gap closing, gap opening, band 
crossing, and band broadening is a generic effect of the 
field modulation. 

When (q, 7) = (2,3), there are eighteen subbands. 
For small values of (3, due to direct touching and indi- 
rect overlapping between the inner sixteen subbands, the 
full energy spectrum exhibits a three subband structure. 
However, there occurs gap opening between the inner 
sixteen subbands with increasing (3, resulting a five sub- 
band structure. We also find that the widths of gaps 
broaden with the increase of (3. For (9,7) = (2,4), there 
are sixteen subbands and they compose a single subband 
structure for small values of /3. However, with increasing 
(3, there occurs gap opening and the full energy spectrum 
exhibits a five subband structure. 



When (q, 7) = (3, 3), there are nine subbands and they 
compose a three subband structure for small values of 
(3. However, gap closing occurs with increasing (3, and 
the three subband structure becomes a single subband 
structure. Further increase of (3 makes the full energy 
spectrum consist of seven subbands. For (q,-f) = (3,4), 
we show the (3 dependence of the energy spectrum in 
Fig. 8, where concurrent occurrence of gap opening and 
gap closing can be seen. In this case, there are twenty 
four subbands and they compose a three subband struc- 
ture for small values of (3. However, as (3 increases, there 
occurs band splitting and the full energy spectrum ex- 
hibits a seven subband structure. Further increase of (3 
makes the full energy spectrum consist of three subbands. 
We summarize in Table 2 the number of subbands of the 
full energy spectra for the parameters (q, 7, (3) taken into 
account. 

The symmetry of the energy spectrum under the 2D 
modulated fields is rather complicated than the cases of 
the ID modulated fields. Se remains for both the 2DSM 
and 2DCM fields as can be seen in Figs. 6-8. However, 
Sx and Sy depend sensitively on the system parameters. 
For the parameters considered, we find that Sx and Sy 
remain for the cases of (q, 7) = (2, 2), (2, 4), (3, 2) in the 
2DSM fields and (q, 7) = (2, 4), (3, 3) in the 2DCM fields, 
while they break for the cases of (9,7) = (2,3), (3,4) in 
the 2DSM fields and (9,7) = (2,3), (3,4) in the 2DCM 
fields. In the meanwhile, the energy spectrum in the 
2DSM field with (q, 7) = (3, 3) has a flipped symme- 
try with respect to k x ^ y ) = 0. Here the flipped sym- 
metry means that the energy spectrum in the range of 
< fc^) < it/R x ( y ) is the same as that in the range 
of —tt/R x ( v ) < k x ( y ) < when E is replaced by — E. 
Sd also depends sensitively on the system parameters; it 
holds (breaks) for the values of (q, 7) that make Sx and 
Sy remain (break). 

Before ending up this subsection, we would like to men- 
tion a few comments. The first is that the above results 
on Se(x,y.d) are not from a theoretical analysis but from 
a numerical study. Thus, further study such as the group 
theoretical analysis 23,26 ' 27 is required in order to deepen 
the understanding of the symmetry breaking. The sec- 
ond is that the k x and k y directions taken into account 
in discussing Sx(y,d) are not the high symmetry direc- 
tions of the Hamiltonian under the 2D modulated fields. 
The reason for taking into account these directions, nev- 
ertheless, lies in testing how Sx(y,d) appeared in the 
Hofstadter's spectrum is influenced by the field modu- 
lation. Since symmetry breaking is generally expected if 
the symmetry axis is not properly chosen, our result that 
the breaking of S X (y.d) depends on the system param- 
eters is quite natural. Finally, it may be worthwhile to 
note the result of Ref. 26, where Shi and Szeto considered 
a Bloch electron in the magnetic field 

B=[B a + {-\) m - n Bi]z (19) 

and found that there is no symmetry breaking in the 
energy spectrum and that the fractal structure remains 
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even though the energy spectrum becomes different from 
the Hofstadter's spectrum, which seems to be contrary to 
our results. However, since the directions used in Ref. 26 
are the highly symmetric (k x ± k y ) directions while the 
directions used in this paper are the k x and k y directions, 
it may not be easy to compare both results directly. The 
only thing we can say is that the field modulation given 
by in this paper [i.e., Eqs. (l)-(3)] is more generalized 
than that of Ref. 26 [i.e., Eq. (19)]; specifying the peri- 
odic magnetic field by a magnetic unit cell, the lattice 
with Eq. (19) becomes the simple checkerboard lattice 
[see, for example, Fig. 2(a) of Ref. 23] while the lattice 
with Eqs. (l)-(3) becomes a generalized checkerboard lat- 
tice [see Fig. 2(b) of Ref. 23]. Thus we expect that the 
latter may present more generic effect of the field modu- 
lation than the former. Further theoretical study on the 
comparison between two cases is also required. 



IV. SUMMARY 



In summary, we studied the effects of spatially mod- 
ulated magnetic fields on the Hofstadter's spectrum. In 
order to obtain rather general properties of the effect of 
the field modulation, we took into account four kinds of 
modulated fields and calculated energy eigenvalues with 
varying q, (3, and 7. Occurrence of gap opening, gap clos- 
ing, band crossing, and band broadening was found to 
be characteristic features of the field modulation, which 
leads distinctive energy band structure from the Hofs- 
tadter's spectrum. It should be noted that these charac- 
teristics are very important since the change of the en- 
ergy dispersion may have crucial influence on the trans- 
port properties; it may allow an electron to change its 
orbit, in view of the semiclassical language. There is 
also a possibility of quenching the integer quantum Hall 
behavior and exhibiting normal behavior when the sub- 
gaps are smeared out. Sensitive dependence of the en- 
ergy spectrum on the system parameters implies that it 
may be difficult to detect direct indication of the Hofs- 
tadter's spectrum in experiment since even a tiny change 
of the magnetic field leads much complicated energy band 
structure. We also studied the effect of the field modu- 
lation on the symmetries appeared in the Hofstadter's 
spectrum and illustrated that the symmetries sensitively 
depend on the system parameters. In this work, we have 
paid attention to the energy spectral properties only with 
a rational flux and the isotropic hopping integral. Since 
an introduction of hopping anisotropy is known to lead 
interesting phenomena like band broadening, 28 it may be 
interesting to study the effect of hopping anisotropy on 
the energy spectrum of a 2D Bloch electron under the 
modulated magnetic fields. 
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FIGURE CAPTIONS 

Fig. 1 Energy dispersion in the uniform magnetic field 
with (a) q = 2 and (b) q = 3. The horizontal plane is 
drawn in units of n. 

Fig. 2 Plot of E versus k y in the 1DSM field with (3 = 0.3, 
where (a) (q,j) = (2,3) and (b) (q,j) = (2,4). The hor- 
izontal axis is drawn in units of n. 



Fig. 3 Plot of E versus k y in the 1DSM field with 
(g,7) = (3,2), where (a) (3 = 0.3, (b) (3 = 0.6, and 
(c) (3 = 0.9. 

Fig. 4 Plot of E versus k y in the 1DSM field with 
(<?,7) = (3,3), where (a) (3 = 0.3 and (b) (3 = 0.9. 

Fig. 5 Plot of E versus k y in the 1DCM field with 
(g,7) = (2,3), where (a) (3 = 0.3, (b) (3 = 0.6, and 
(c) (3 = 0.9. 

Fig. 6 Energy dispersion in the 2DSM field with (q, 7) = 
(2, 2), where '(a) (3 = 0.3 and (b) (3 = 0.9. 

Fig. 7 Plot of E versus k y in the 2DSM field with 
(q,j) = (3,2), where (a) (3 = 0.3, (b) (3 = 0.6, and 
(c) (3 = 0.9. 

Fig. 8 Plot of E versus k y in the 2DCM field with 
{q, 7) = (3, 4), where (a) (3 = 0.3 and (b) (3 = 0.6. 



TABLE CAPTIONS 

Table 1 Number of subbands of the full energy spectrum 
under the ID modulated fields. 

Table 2 Number of subbands of the full energy spectrum 
under the 2D modulated fields. 
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